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Abstract 



The paper presents a generalization of Arnold-Falk-Winther elements for three dimensional 
linear elasticity, to meshes with elements of variable order. The generalization is straightforward 
but the stability analysis involves a non-trivial modification of involved interpolation operators. 
The analysis addresses only the /i-convergence. 
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1 Introduction 

Linear elasticity is a classical subject, and it has been studied for a long time. The paper focuses on the 
so-called dual-mixed formulation with weakly imposed symmetry that may be derived by considering 
stationary points of the generalized Hellinger-Reissner functional lfl6l . We restrict ourselves to the 
static case only and, for the sake of simplicity, we assume that the body is fixed on the whole boundary. 
We look for stress tensor a G H(div, 0;M), displacement vector u G L 2 (Q; V), and infinitesimal 
rotation p G L 2 (f2;I€) satisfying 



The first equation represents a relaxed form of the Hooke's law combined with Cauchy geometrical 
relations, the second one represents the equilibrium equations (in a strong form), and the third one 
enforces the symmetry of the stress tensor. We refer to the next section for a detailed description of 
energy spaces: H(div, M), L 2 (0; V) and L 2 (f2; K). The operator A denotes the compliance tensor 
(operator) mapping stress tensor into strain tensor. The operator is bounded, symmetric, uniformly 
positive definite, and it preserves the symmetry of the tensor. 

The traditional motivation for studying the formulation (11.11) comes from handling nearly incom- 
pressible materials. Our interest in the subject stems from a study of a class of (visco)elastic vibration 
problems for structures with large material contrast, see ifTTl for a motivating example. 

A number of authors have developed approximation schemes based on formulation (11.11 ). among 
others seeSI21Ii21IIlEll3II31l3IISllISllEll2aED. For a brief description of these methods, we 
refer to the introduction in |5]. We also refer to the recent work of Cockburn, Gopalakrishnan and 
Guzman [ 8 ] who have developed a new mixed method for linear elasticity using a hybridized version 
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The work presented in this paper is based on the mixed finite element methods developed by 
Arnold, Falk and Winther in Bl[T2ll5l. 

The ultimate goal of this work is to lay down theoretical foundations for, and implement a fully 
automatic fop-adaptive Finite Element (FE) method based on a generalization of the AFW element to 
meshes with variable order. The generalization builds on the exact grad-curl-div sequence that holds 
for hp meshes, see ifTUl [TT il and it is rather straightforward. The formulation is easily accommodated 
in a general hp code supporting the exact sequence. 

The h convergence analysis presented in Q for meshes with arbitrary but uniform polynomial 
order, does not however generalize immediately to elements with variable order. 

With the proof of p and hp convergence as an ultimate goal, our initial efforts start with a less 
ambitious goal of proving first stability and convergence for uniform /i-refinements of meshes of 
variable order. 

At the first glance, a generalization of the techniques from [5 ] seems to be easy. But, as we have 
shown in ifTTl , the (natural generalization of) canonical projection operators defined in |@1 do not 
commute with divergence operator on meshes with variable order, a property essential in the proof of 
discrete stability. We have resolved this problem by invoking the Projection Based (PB) interpolation 
operators in IfTTl . Unfortunately, the PB operators do not commute with an algebraic operator S n -2, 
introduced in another essential construction in the AFW proof of discrete stability. We have 
resolved the problem by designing a new, special operator Wh in IfTTl that satisfies the commutativity 
property, as needed. Unfortunately, we managed to prove well-definedness of operator Wh only for 
polynomial orders < p < 3, and only for two space dimensions, see IfTTl . In this contribution, 
we resolve the problem by designing new PB operators and a a variant of Wh, a new operator H f \ , 
discussed in the text. 

An outline of the paper is the follows. Section 2 introduces notations. In Section 3, we define 
the involved finite element spaces. Section 4 recalls the mixed formulation of linear elasticity with 
weakly imposed symmetry and states the Brezzi conditions for the stability. In Section 5, we establish 
all technical results needed for proving the stability. We construct the new PB operators and the 
operator Hf\ . Finally, in Section 6, we prove the Brezzi conditions. 

2 Notations 

In this section, we introduce some basic notations. We define M to be the space of 3 x 3 real matrices, 
and V to be R 3 . For any 3x3 real matrices A, B, we define 

A: B = tr(AB T ). 

We denote by S and K the subspaces of symmetric and anti-symmetric matrices in R 3x3 . Each anti- 
symmetric matrix can be identified with a vector in V given by the mapping vec : K — > V: 
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$7 is a domain in R 3 . For any vector space X with inner product, we denote by L 2 (£l; X) the space 
of square-integrable vector fields on Q with values in X. In the paper, X will be R, V, M, or IK. 
When X = R, we will write L 2 (Q,). The norm associated with L 2 (Q,; X), denoted by || • Hz^^.x), is 
obtained by taking the square root of the sum of (squared) I? norms of individual components of the 
vector fields on Q,. 
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Notice that, for scalar- valued functions, norm || ■ ||.L 2 (n) coincides with the standard L 2 -norm. The 
corresponding Sobolev space of order m, which is the subspace of L 2 (Sl; X) consisting of functions 
with all partial derivatives of order less than or equal to m in L 2 (Sl; X), is denoted by H m (Sl; X). 
The norm associated with H m (S2; X), denoted by || • \\H m (tt-,x.)> equals to the square root of the sum of 
(squared) L 2 -norms of all partial derivatives with order less than or equal to m, for all components of 
vector fields on 11 When X = R, || • coincides with the standard i/ m -norm for scalar-valued 

functions. 

The spaces .ff (curl, O), H(div, SI) are defined by 

# (curl, SI) = {u£ L 2 (Sl;Y) : curki G L 2 (Sl;Y)} 
H(div,n) = {ve L 2 (S1;Y) : divv G L 2 (S1)} 

with the norms, 

\\u\\h(cw\,Q) = (\\ u \\h(n-Y) + l|curlw||| 2(n . v) ) 1/2 , u G #(curl, SX) 

||w||H(diy,n) = (IMl£a(n ; v) + \\ d ™ v \\h(n)) 1/2 i v G H (tot, SI). 

We extend the definitions of V to V-valued functions, curl and div to M-valued functions by 
applying these operators row-wise. The space //(curl, S}; M), if (div, SI; M) are defined by 

fl"(curl,fi;M) = {a G L 2 (Sl;M) : curia G L 2 (Sl;M)} 
F(div,f);M) = W € L 2 (S1;M) : diva G L 2 (S1;Y)} 

with the norms, 

lklltf(curl,J2;M) = (IMIl^M) + I l curlcr 1 1 is (s^ ; MI) ) 
IkllH(divAM) = (lkllL2(n;M) + ll divcr llL2(n ; V)) 1/2 - 

V r (ST) denotes the space of polynomials on SI with degree less than or equal to r. When r is a 
negative integer, V r (ST) = {0}. V r (Sl; V) = [P r {Sl)} 3 . Throughout this paper, we assume that r is a 
nonnegative integer. 

3 Finite element spaces 

3.1 Finite element spaces on a single tetrahedron 

Let T be an arbitrary tetrahedron in R 3 . We denote by Afc(T), the union of /c-dimensional subsim- 
plexes of T. We denote by A(T), the union of all subsimplexes of T. 
For any r G Z + := {n G Z : n > 0}, we introduce 

V r A 3 (T) := V r (T),V r A 2 (T) := V r (T;Y), 

V r A 2 (T) := {u G T r A 2 {T) : VF G A 2 (T),w • n\ F = 0}, 

V~A 2 (T) := Vr-i(T;Y)+xV r -i(T), 
V;A 2 (T) := {u G P~A 2 (T) : VF G A 2 (T),u • n\ F = 0}, 
P-A^T) := V r -i(T;Y) + x x V r -i(T;Y), 
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l>~k l {T) := {co G V~K l {T) : VF G A 2 (T),u - (w • n)n| F = 0}. 

V r A°(T) := {u G P P (T) : u| 9T = 0}. (3.1) 

Here n is a normal unit vector on F. For F, an arbitrary face of T, we introduce 

V-A^F) := Vr-iiF-^+yVr-^F). (3.2) 

Here y denote any orthogonal coordinates on F. In HOG!, spaces in (13.113.21 ) are defined in the 
language of exterior calculus. Here we rewrite them in the standard language of calculus. Please refer 
to El and HP for a detailed correspondence. 

We denote by f a mapping from A(T) to Z + such that if e, / G A(T) and e C / then f(e) < 
We introduce now formally the FE spaces of variable order. 

Definition 1 

V f A 3 (T) := P f{T) A 3 (T) = V H T)(T), 
VrA 2 (T) := {u G V r , {T) A 2 (T) : VF G A 2 (T), w • n| F G P r - (F) (F)}, 
PrA 2 (T) := G Pr (T) A 2 (r) : Vi? G A 2( r ). w • n| F G ^ (F) _ 1 (F)}, 
PrA^T) := {u G n"(T) Al ( T ) : Vi? G A 2(n w - (u • n)n| F G P^A^F); 
Vt G Ai(T), w • t G Pf( e )_i(e)} where t is a tangential vector on e. 

Remark 2 In the definition of A l {T), u — (to ■ n)n\p is a tangential vector field on F. So 
oj — (oj ■ n)n\p can be considered as a two-component vector field. 

Definition 3 We define V f A 3 (T;V) := V f ( T ) (T; V). We also define VfA 2 (T; V),VfA 2 (T; V), and 
VJ A 1 (T; V) as matrix-valued polynomial spaces whose rows stay in Vf A 2 (T), V~A 2 (T), andV-A 1 ^) 
respectively. 

Remark 4 Finite element spaces in definitions\l\and\3\are the same as those introduced in [17] for 
n = 3. In this paper, we just rewrite them in the language of standard calculus. 

3.2 Finite element spaces on a bounded polyhedral domain 

Let Th be a tetrahedral mesh. Here h represents the biggest diameter of tetrahedrons in Th- We extend 
the map r to a mapping from A (7^) to Z + such that if e C /, then f(e) < r(f). For any T G A^(Th), 
the restriction of r to A(T) is represented with the same symbol f. We denote by A^{Th) the union 
of fc-dimensional subsimplexes of Th, and by A(Th) the union of all subsimplexes of Th- 

Definitions Let T h be a tetrahedronal mesh. We define VfA 3 (T h ), VfA 2 (T h ), V~A 2 {T h ), V~A 1 (T h ), 
VfA 3 (%;Y), VfA 2 (Th;Y),VfA 2 (Th;Y), and VfA 1 (T h )Y) as spaces of piece-wisely smooth func- 
tions or vector fields on Th whose restrictions on T are VfA 3 {T), VfA 2 (T), VfA 2 (T), V~~A (T), 
P r -A 3 (T;V), VfA 2 {T;V),V~ A 2 (T;V), andV~A x (T; V) respectively, for any T G A 3 (7h). 

Remark 6 Obviously, we have 

VrA 2 (Th),V~A 2 {Th) C H{div, n), V^A\Th) C H(curl, O), 

V f A 2 {Th\V),V~ A 2 {T h ]V) C H(<fiv,n;M), VjA 1 {T h ;N) C ff(c«r/,0;M). 

//ere f2 ij aw o/?e?2 subset in R 3 w/f/z = Ure^F. Tfte spaces defined in (T5]) /jave fteen introduced in 
HI 71/ using the language of differential forms. 
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4 Algebraic operators and some auxiliary properties 

In this section, we will introduce two algebraic operators, and prove some of their relevant properties. 
Definition 7 We introduce a linear map S2 defined as follows, 

S 2 U = (u 23 - U32, u 3 i - U13, U12 - u 2 i) T ■ 
Here U is an arbitrary matrix in M 3x3 . 
Remark 8 It is easy to check that S2U = vec(U T — U). 

Definition 9 We define linear map S\ as follows, 

S X W = W T - tr(W)I. 
Here W is an arbitrary matrix in R 3x3 . 



Lemma 10 Operator Si is invertible. And S 1 1 W = W T — ^tr(W)I. 

Lemma 11 divS\W + S 2 curlW = 0, MW G H X {Q, M). Here, Q is any open subset in R 3 . 
Proofs of Lemma [TOl and Lemma[TT]are straightforward. 

Lemma 12 Let Tbea tetrahedron in M 3 , and W £ H l (T, M). Let F be any face ofT. If W-t\ F = 
for all tangential vectors t on F, then S\W ■ n\p = where n is an unit normal vector on F. 

Proof. According to definition |9j we have 

-W22 ~ W 33 W 2 1 U> 3 1 

S\W ■ n = W12 -w\\ - u>33 u>32 • n 

W13 W23 —Wn — W22 

= [nx(-w 2 2 ~ W33) + n 2 W2i + 7i 3 w 31 , 7111012 + n 2 (-wu ~ w 33 ) + n 3 w 32 , 
niw 13 + n 2 w 2 3 + na{-w\\ - W22)] 1 ■ 

In the following, we will show that n\{— W22 — W33) + n2W2\ + ^.3^31 = on F . The proof of 
the other two components of S\W ■ n being zero on F is similar. 

Obviously, (n 2 , — n\, 0) T -n = and (713, 0, — ni) T -n = 0. This implies that W-(n 2 , —n\, 0) T = 
and W ■ (ns, 0, — rii) T = on F. So we have —711W22 + n 2 W2i = and —711^33 + re 3^3i = on 
F. This shows that n\{— W22 — W33) + 71-2^21 + ^3^31 = on F. m 



Lemma 13 For any W, Q 6 M, we have S\W : Q = W : S\Q. 
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Proof. 



S X W : Q = 



-W 2 2 ~ W 33 W 2 l 

Wu -Wu - W 33 



Wl 3 



W 23 



W 3 l 
W 32 
-Wn - W 2 2 



=(-W22 - w 33 )qn + W2iqu + w 3 iqi 3 

+ ^12921 + (-Wll - W 33 )q 2 2 + ^32923 
+ ^13931 + 1^23932 + {-Wn - W22)<?33 
=w n (-q 2 2 ~ q 33 ) + ^12921 + ^13931 
+ ^21<?12 + W 2 2(-qil ~ q 33 ) + 1^23932 
+ W 31 q 13 + 1^32923 + W 33 (-q n - 922) 

— ?22 — 933 921 931 

- : 912 -911 - 933 932 

913 923 _ 9ll — 922 



Q 



W : 5iQ. 



Lemma 14 Let T be a tetrahedron in M 3 . We take W G P^A^T; v )- V 

Js 1 W:Q = 0, Q € V r -i(T;M), 



then W = 0onT. 

Proof. According to Lemma [T3l we have 



J W:S 1 Q = 0, QeP M (T;l 



(4.1) 



By the definition of Si, it is easy to see that SiV r -i(T; M) C V r -i(T; M). According to Lemma 
[TOl we conclude that SiV r -i(T; M) = V r -i(T; M). According to Lemma 4.11 in [4], we have that 
W = 0onT. m 



5 Mixed formulation for the elasticity equations with weakly imposed 
symmetry 

We begin by rewriting formulation (11.11 ) using operator S2. The elasticity problem becomes: Find 
(a,u,p) G H(div,n-M) x L 2 (ft;V) x L 2 (S1;V) such that 

(Ao-,t) + (divr,?/) - (S 2 T,p) = 0, t G -ff (div, Q; M), (5.1) 
(<fiv<7,«> = (/,«), t)a 2 (fi;V), 
(S 2 a,<z} = 0, gGL 2 (Q;Y). 

Here (■, •) is the standard L 2 inner product on fl This problem is well-posed in the sense that, for 
each / G L 2 (fi; V), there exists a unique solution (a,u,p) G ff(div,fi;M) x L 2 (S1;V) x L 2 (S1;V), 
and the solution operator is a bounded operator 

L 2 (fi;Y) — ► #(div,f};M) x L 2 (ft;V) x L 2 (S1;V). 
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See El and flU for the proof. 

Next, we consider a finite element discretization of (15.11) . For this, we choose families of finite- 
dimensional subspaces 

A 2 h {M) C #(div,ft;M),A^(V) C L 2 (fl; V), Aj(V) C £ 2 (^;V), 

indexed by /i, and seek the discrete solution (<Jh,Uh,Ph) E A|(M) x A|(V) x A h (Y) such that 

(Ar ft ,r) + (divr I u fc )-(52r,p fc > = 1 r G A|(M), (5.2) 
<div(7 ft ,t;> = (/,«>, «eA£(V), 
(S 2 a h ,q) = 0, 

The stability of (15.21 ) will be ensured by the Brezzi stability conditions: 

(SI) ll r llH(div,Q;M) ^ c i(^ r > r ) whenever r G A|(M) satisfies (divr,u) = (5.3) 
Vi; € A| (V) and (S 2 r, g) = Vg G Aj(V), 

(S2) for all nonzero (v, q) G A^(V) x A /l (V), there exists nonzero (5.4) 

r G A|(M) with (divr,v) - (S 2 r,q) > c 2 ||T||^ (divAM) (|H| i 2 (n . v) + ||g|| £ 2 (n . V )) > 

where constants c\ and c 2 are independent of h. 

For meshes of arbitrary but uniform order, conditions (15.31 ) and (15.41) have been proved in Q 
and [12]. In what follows, we will demonstrate that they are also satisfied for meshes with elements 
of variable (but limited) order. In this paper, we define A|(M) = Vf+\A 2 (Th]Y), and A|(V) = 
A^(V) = VfA 3 (%;'y). We assume that there is r max G N such that for any h > and / G 

I max- 

6 Preliminaries for the proof of stability 

From now on, we assume that is a bounded polyhedral domain in M 3 . We also use the standard 
assumptions for shape regular meshes, which means that the ratio between outer diameter and inner 
diameter of any tetrahedron in any mesh has an uniform upper bound. 

In the proof of stability, the following three commuting diagrams are essential. 



uj h | n? h | (6.i) 

Vr +1 A 2 (T h ;V) ^ VrA 3 (T h ;V) 
H^M) ^ L 2 (n ; Y) 

(6.2) 

Vi +1 A 2 (T h ;V) "H^ VrA 3 (T h ;V) 

I n^;~ | (6.3) 
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Here IT- h is the I? orthogonal projection operator onto VfA 3 (Th,;Y). II? h , n~' h , and H"-^ are 

projection operators into Vf+\A 2 (Th]Y), P f+1 A 2 (7/ i ; V), and V^ +2 A l (Th', V) respectively. 

In |H|5l, the canonical projection operators introduced by Arnold, Falk and Winther can make 
(I6.1I6.21I63T ) commute for meshes with uniform order. But for meshes with variable order, the natural 
generalization of the projection operators fails to make commute both d6.ll ) and ( 16.21 ), see a counter- 
example presented in the appendix of ifTTl . To overcome the difficulty, we recalled Projection Based 
(PB) interpolation operators from IfTTl . According to Lemma 19 and Lemma 20 in IfTTl . there exists 
projection based interpolation operator II? h , which satisfies the following properties. 

dMlf ifc T = n? h divr, r G H\Sl;M). (6.4) 

\\Ul h r\\<C\\T\\ m{Q . u) , reff^flii). (6.5) 

Here C is independent of r and h. Though Lemma 20 in IfTTl has been proved for quasi-uniform 
meshes only, it is straightforward to extend it to get ( 16.51 ) for shape regular meshes as well. Please 
refer to ifTTIl for the details on the PB interpolation operators. 

With the PB operators in place, the difficulty shifted to to defining a special projection operator 
Hfh , denoted by Wh in IfTTl , that makes now (16.3b commute. The commutativity property followed 
directly from the construction of Wh but proving that it is well-defined, turned out to be difficult. We 
managed to show only that is well-defined for < r < 3 with n = 2. In the following, we 

2 — — 1 — 

will use a different reasoning to demonstrate that there exist projection operators ILA and Iiy^ , both 
well-defined, that make (16.2U6.3I) commute for arbitrary 3D meshes of arbitrary order. Note that the 
operators will not be constructed explicitly. 

6.1 Projection operators on a reference tetrahedron 

Let T be a fixed tetrahedron in M. 3 . We are going to design projection operators H^ and IT _'T into 
V~ +l k 2 {f- V), and V~ +2 h l {f- V) respectively. 

Definition 15 We take r to be a mapping from A(T) to 7L+ such that ife,f € A(T) and e C f, then 
r(e) < f(f). We put k = dim curl & V~,f )+1 A 1 (T;V). 

We define {/ P>1 , • ■ • , / f fc } as a basis of cwrZ^-^^A^T; V). We define {g fA , ■ ■ ■ ,cfr,k} 

as a linearly independent subset ofV~^ _i(T; M) such that {fff^i, • ■ ■ , 9?,k} © ^£^ > rA 3 (T; V) = 

P f » f ,. 1 (f;M). We define h f>i (x,t) = (1 - t)f f {x)+tg r ~.{x)forany 1 < % < k. 

Remark 16 It is easy to check that k = dimP p(:f , ) _ 1 (f;M) - dim V i 7 ? fA 3 (T ; V). Take r = r(f). 

Then dim curl & V r +i A 1 (T;Y) = 3(dimP r+ iA 1 (f ) - dim£ r+2 A°(f )) = ^(2r + 5)r(r - 1). 

dimP r _i(f ;M)-dimV ii; P r A 3 (f ;V) = 3(dimV r -i(f ■ V) - dim P r (f )/M) = -(2r+5)r(r-l). 
For f/je dimensions of finite element spaces mentioned above, please refer to formula (3.1) in Ml/ aid 
page 51 m Ml/- 
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2 i / * 

Definition 17 For any t G [0,1], we define the linear operator TlJ^^ mapping H (T;M) onfo 
"P^ 1 A 2 (T 1 ; V) by the following conditions. 

fdiv & (U^ t U-U)-fjdx = 0, »)Gp f(f) (f;V)/M. (6.6) 
/ (n 2 'T U-U): hfiix. t)dx = 0, 1 < i < k. (6.7) 

Jj, r,T,t ' N ' 

Jjnl'^ t U-U).h}-jlds = 0, FgA 2 (T), deV f{P) (F;V). (6.8) 

Definition 18 For any t € [0, 1], we define the linear operator nl'T ^ mapping H 1 (T; M) into 
'P^ 2 A 1 (T; V) Zry the following conditions. 

[div & S 1 (nl'ZjV-W)-fjdx = Q, T}€V f( fi(f;V)/R. (6.9) 

/ 5i(nl'T W - W) : f)da; = 0, 1 < < k. (6.10) 

jMl'^W-W)-i}-fids = 0, F£A 2 (f), /i6^,(F;V). (6.11) 



nl'+W-t\, = 0, eGAi(T). (6.12) 



f,T,t 



In (16.81 ), n is a unit normal vector on F. In ( 16.111 ), i is any tangential vector on F. Notice that the 
dimension of tangential vector space on F is two. In (16.12I ). t is a tangential vector along e. 



Lemma 19 For any r(T) G Z+, operator n 2 '^, ^ jj a linear projection, and a well-defined operator 
for all but finitely many values oft G [0, 1]. 

Proof. It is easy to see that the conditions (16.616.716.81 ) are well-defined for any U G H 1 (T;M). 

2 — 

Obviously, if n_'^ ^ is well-defined, then it is linear and a projection. It is sufficient to show that for 

any U G V~ +l k 2 {f; V),U = if tf'T U = 0. 

By Theorem 4.12 in £Q, U G V^ f)+l A 2 (f;Y) because 

f[U-n]-M$ = 0, FeA 2 (T), fiEV- {P} {F;W). 

J F 

So it is sufficient to show that, for any U G "P^^A^T; V), U = 0, provided, 

J^div^U ■ f/dx = 0, j)GP f(f) (T;V)/l; (6.13) 
J U:hf ji (x,t)dx = 0, 1 <i < k. (6.14) 
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Since U G V f ^ +1 A 2 (f; V), (16.131 ) can be integrated by parts to yield 

Jjj:V S: fjdx = 0, fjeP f{f) (f;V). (6.15) 

According to the definition of hf^x, t) and the fact that "P r ~A 3 (T; V) = ~Pfrf\ (T; V), the assertion 
is true for t = 0. Indeed, when t = 0, conditions (16.1516.141 ) can be rewritten as 

j^U:Qdx = 0, geP^ff;!), 
By Lemma 4.11 in |@), we have £/ = 0. This implies that n 2 'T is well-defined for t = 0. 



We denote by C(i) the matrix associated with the left hand side of conditions (16.616-716- 8b - Then 
n 2 'T is well-defined if and only if C(t) is a square nonsingular matrix. We have already known 

that C(0) is a square nonsingular matrix. So C(t) is a square matrix for any t G [0, 1]. Notice that 
det(C(i)) is a polynomial of a single variable t. Since det(C(0)) ^ 0, then there are at most finitely 
many t G [0, 1] which make det(C(i)) = 0. This implies that n 2 'T is well-defined for all but finitely 
many values of t G [0, 1]. ■ 

Lemma 20 For any f(T) G Z+, operator IIjT ^ jj a well-defined, linear projection operator for all 
but finitely many values oft G [0, 1]. 

Proof. It is easy to see that the conditions (16.916. 1016. Ill 16.121) are well-defined for any € 
i/ 1 (T;M). Obviously, if EL'T is well-defined, then it is linear and a projection. It is sufficient 

to show that for any W G T^A^T; v ) witn W" • t|e = 0, then W = if nl'T TV = 0. Here e is 

any edge of T, and i is a tangential vector along e. 

By Theorem 4.12 inft^G V~ ^ A^f ; V) because 

f[W-t]-fids = 0, FgA 2 (T), fi€V f ^(F;V); 

J F 

W-i\e = 0, eeAi(f). 
So it is sufficient to show that, for any W" G ^ja +2 -^ 1 (^'! V )' ^ = °> P rov i ded > 

J div & S 1 W-f)dx = 0, i) eV Hf) (f;Y)/R; (6.16) 

/ SiW" : hf,i(x,t)dx = 0, 1 < i < Jfe. (6.17) 

Notice that W" G PI^^A^T; V). By Lemma [TU we have that 5j W" ■ n|^, = for any F G 

A2(T). So we can integrate thus (16.161 ) by parts without obtaining any boundary term. Condition 
(16.161 ) can be rewritten as follows. 

Js 1 W:V & fjdx = Q, f)£V fif) (f;Y). (6.18) 
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By the definition of hp t i(x, t), for t = 1, conditions (16.181) and (16.171) can be rewritten as 

JSxW :Qdx = 0, QeP^ H (f;M). 

Lemma [141 implies then that W = 0. This shows that IIjT ^ is well-defined for i = 1. 

We denote by C(i) the matrix associated with the left hand side of conditions (16.916. 1016. 11 16. 121 . 
Then EL'T ^ is well-defined if and only if C (t) is a square non-singular matrix. We have already known 
that C(l) is a square non-singular matrix. So C(t) is a square matrix for any t G [0, 1]. Notice that 
det(C(i)) is a polynomial of a single variable t. Since det(C(l)) 7^ 0, then there are at most finitely 
many t G [0, 1] which make det(C(i)) = 0. This implies that ^ is well-defined for all but finitely 

many values of t G [0, 1]. ■ 

According to Lemma [191 and Lemma l20l we can choose t r G [0, 1] for any f(T) such that both 
IT'I and IT'I are well-defined. Here t r depends only on r(T). 

2 2 1— 1 

Definition 21 We <i<?/»i<? operators H ' ~ := IT ' . arcd IT ' - := II ' - . 

f,T r,T,t r r,T r,T,t r 

6.2 Projection operators on a physical tetrahedron 

Let T be an arbitrary tetrahedron in R 3 . Then there exists an affine mapping from the reference 
tetrahedron T to T, defined by 

x = Ax + b. (6.19) 

Here A is a 3 x 3 real non-singular matrix, and b is a vector in M 3 . In the following, we always relate 
x and x by (I6.19I ). We take f to be a mapping from A(T) to Z + such that if e, / G A(T) and e C /, 
then f(e) < r(/). In the following, we denote by x(x) the inverse of the affine mapping described 
above. 

Definition 22 We define the linear operator II?' T mapping H l (T]~M) onto Vf +l h?{T;V) by the 
following conditions. 

/ div(nf~U-U)-7]dx = 0, j) G Tft T) (T;Y)/R. (6.20) 

JT 

J (U 2 -'~U - U)(x) : [Ahf }i (x(x),t r )A~ 1 ]dx = 0, 1 < i < k. (6.21) 

[ [(Ul'-U - U) ■ n] ■ pds = 0, FGA 2 (T), p eV f(F) (F;Y). (6.22) 

Definition 23 We define the linear operator ni'^ mapping H 1 ^; M) mto 7 3 ^|_ 2 ^ 1 (^> ^0 by the fol- 
lowing conditions. 

I divSi(Ill'~W - W) ■ rjdx = 0, 77 G V f{T) (T; V)/M. (6.23) 
/ SiCnJ'yW-W) : [Afe ?)i (a;(a;),i r .)A" 1 ]da; = 0, 1 < i < k. (6.24) 
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[ [{nl'~W - W) ■ t] ■ fids = 0, FgA 2 (T), fjL€V m (F)\). (6.25) 
Jf 

ul'-W-t\ e = 0, eeAi(T). (6.26) 

We want to "pull back" f from A(T) to A(f ). We put f(e) = f(e) for any e G A(f ). Here 
e := Ae + b. Then we have the following lemma. 

Lemma 24 For any U,W G fl^T; M )> w<? U,W £ H x {f; M) fry 

i7(as) = A" T ?7(a;)A T , PF(a;) = AW^-A -1 . 

Then we have 

U 2 f '-U(x) = A- T Il 2 _'fU(x)A T , ul'~W(x) = AUl^W(x)A-\ xeT. (6.27) 

2 — 1 — 

So operators n~' T and n~' T are well-defined. 

Proof. For the result of n~' T , the proof is straightforward. For the result of Hf' T , we need utilize the 
definition of Si. Notice that 

SiW{x) = W{x) T - ti(W(x))I = A- T [W{x) T - ti(W(x))I]A T = A~ T SiW(x)A T . (6.28) 

Using (16.281) . it is now straightforward to prove the result for IT-'^ . ■ 

Lemma 25 For any U G H l {T\ M), we have 

Uf T divU 2 J~U = H$ T divU. 
Here II? T is the orthogonal projection operator from L 2 (T; V) onto VfA 3 (T; V). 

Proof. According to the definition of 11-^ , we have (I - n?^)n^C7 = for any U G H X {T\ M). 
So it is sufficient to show that IT? T divU = for any U G H 1 (T; M) with ILf^U = 0. 

Now, we choose U G H X (T; M) with U^U = 0. We only need to show that j T divll ■ fj = 
for any rj G V?(T)(T;'V). Obviously, we can choose c G M 3 such that rj = rj + c, where rj G 
P f ( T ) (T; V) /R. Then we have 

/ divll ■ rjdx = / divll ■ r/dx + / divll ■ cdx = / divll ■ r/dx + / (U ■ n) ■ cds. 
Jt Jt Jt Jt Jot 

By (16.20l ). (16.22l ) and the face that 11-'^ J7 = 0, we have f T divU ■ fjdx = 0. This implies that 
n| T divn?'~C7 = n? T divC/ for any U G H 1 ^; M). ■ 

Lemma 26 For any W G H X (T\ M), we have 

U^-SiUl'-W = uf^SiW. 
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Proof. According to the definition of n~' T , we have 

(J - n£~)n~;~ w = o, weh\T;M). 

So it is sufficient to show that S{W = for any W G F X (T; M) with = 0. 

Now, we choose W G F 1 ^; M) with Il£~W = 0. By (16331) and (f6T24l >. we have, 

J divSiW -rjdx = 0, 7] G P f(T) (T; V)/R; 

5iW : [^^(scjt,.)^ -1 ]^ = 0, 1 
In order to demonstrate that H~ T S\W = 0, we only need to show that 

[SiW-n] -fids = 0, FGA 2 (T), /iG? f(f) (F;V). 



X 



According to the definition of Si, we have, 
S{W-n = 



-W22 ~ W 33 
W12 
W\ 3 



1021 
-Wn - W 33 

^23 



1031 
^32 
-Wu - W 2 2 



n 



(n 2 W2i - niw 2 2) + (n3W3i - niw 33 ) 
-(n 2 w n - mwi 2 ) + (n 3 w 3 2 - n 2 w 33 ) 
-(n 3 w n - niw 13 ) - (n 3 w 2 2 - n 2 w 23 ) 

Consequently, for any ji G V?tp\{F; V), we have, 



[S 1 W-n]-fi= W 





U2 




" -1*2 ' 




-ni 


) 


Ml 















n 3 




" — A*3 " 

















) 





+ i^w ■ 


™3 


) 


-Ma 








Mi 




_ -n 2 _ 




M2 



By (16.251 ) and the fact that nl' T PF = 0, we conclude that 

f [SxW ■ n] ■ ixds = 0, FGA 2 (T), /xGP r ~ (F) (F;V). 
Consequently, n^ySiW = 0. ■ 

Lemma 27 77jere exists c > smc/j ?/jatf, /or any U, W G .ff^T; M), 

l|n r V T C/||L2( T;M ) < c||f/|| H i( T;M ); (6.29) 

||cMr/n~;~VF|[ L 2 (T;V) < c(/i^ 1 1| W\\ L 2 {T . M) + \\W\\ H i {T . M) ). (6.30) 

//ere hx is the outer diameter of T, and c is independent of T. The constant c may depend upon the 
ratio of outer and inner diameters ofT. 

Proof. (16.29l6.30l are obtained by standard scaling techniques. The proof for ( 16.291 ) is the same as 
that for Lemma 20 in ifTTl . The proof for ( 16.301 ) is the same as that for Lemma 29 in IfTTl . ■ 
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6.3 Projection operators on tetrahedral meshes 

As we stated at the beginning of this section, we use standard assumptions for regular meshes. This 
means that the ratio between outer diameter and inner diameter of any tetrahedron in any mesh has a 
uniform upper bound. We are going to extend operators n~' T and n~' T now to the whole mesh Th in 
such a way that they make ( 16.216.31 ) commute. 

Definition 28 We define mappings IT?'" : H x (Sl;M) Vf +1 A 2 (T h ;^) and ul \ : ff^fyM) -> 
V~ +2 A\T h ;V)by 

(n r 2 ; -f/)| r = n?;-(c/| T ) ; (nJ;-^)| T = nl>-(w\ T ). 

Here T G A 3 (7^), and 17, W G M). 



Lemma 29 For any U,W £ H^OjM), n~;~Z7 € 7y +1 A 2 (Th;V) and € P^A 1 ^; V). 

And we /jave 

i%divn?£[/ = ul h divU; u^s^w = u^s.w. (6.31) 

And there exists a constant c > 0, which is independent ofTh,U,W, so that 

\\^f , ,h U \\L 2 (n;M) - C W\\H 1 (n;M)'i (6.32) 

||cwr/n~^W| T || L2(T;V) < c(/i^ 1 ||Ty|r||i2 (T;M) + ||W|r||/ri(T;M))- ( 6 - 33 ) 
Here T G A3 (7^), and Iit is the outer diameter ofT. 

Proof. This is by definitions of n?'^ and n~'^, Lemma l25l Lemma l26l and Lemma l27l ■ 

To remove the Kf} factor in ( 16.331 ). we introduce a Clement-type interpolant Rh mapping H 1 (£1; M) 
into continuous piece- wise linear M-valued function on Th (The operator H h in Theorem 5.1 of Q, 
using example 1 in Q). Then there exists a constant c > such that \/W G M), T G Th, we 

have that 

\\W - RhW\\ L 2(T-M) < C^tII^Hhi^M); \\ W - R hW\\ H i(T;M) < c||W|| fl -i( ST . M ). 
Here := \Jt' eT h -T' nT^q> T ■ Then we follow [6] and define Ii^ h := Up h (I - R h ) + R h . 
Lemma 30 n~£ maps from H 1 (Q; M) into Vf^A 1 {%) V). 

n 2 ;-5in^w = uf-s 1 w, weH\n;M). (6.34) 

And there exists a constant c > .rac/z that for any W G il (f2; M), 

llcwrZni'^ VF|| L 2 (n . v) < cllWHjjifQ-M). (6.35) 

Proof. Since Rh maps if 1 (r2;M) into continuous piece-wise linear M-valued function on Th, we 
have n~'^ maps from H 1 (Q; M) into V^A 1 {%] V). The proof for (16.3416.351 ) is straightforward. ■ 
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7 Stability of the finite element discretization 

We will use the following well-known result from partial differential equations, see Ifl4l . 

Lemma 31 Let Q be a bounded domain in IR 3 with a Lipschitz boundary. Then, for all fj, G L 2 {Q), 
there exists 77 G H (T2;V) satisfying divr] = fi. If in addition, f^pdx = 0, then we can choose 
rj G H\n;Y). 

Remark 32 The domain Q need not be contractible. 

The main result of this paper is the following theorem. In the proof we follow the lines of proof 
of Theorem 9. 1 in [12] , Theorem 7. 1 in Q and Theorem 1 1 .4 in |4| . The main difference is in the use 
of our operator ny ft in place of the operator ft^~ 2 from f!2l . 

Theorem 33 Let Vibe a bounded polyhedral domain in R 3 with a Lipschitz boundary. We assume 
that the meshes are regular. Then for any G VfA. 3 (Th;Y) x VfA 5 (Th]Y), there exists a G 

Vr+i^iTh', V) such that diva = \i, —II 3 j-_ S2G = w. Arcc? we have 

\\c\\H(div,n;M) < c (IMU 2 (n ; V) + IHIl^CjV)), ( 7 -l) 
where the constant c is independent ofuj, fi and h, but it may depend upon max rg ^ 3 ^ ) r(T). 

Proof. We want to show that Brezzi stability conditions (15.3I ). (15.4I ) are satisfied. The condition (15.31 ) 
is obviously satisfied since, by construction, cUVPf + iA 2 (f&; V) C VfA. 3 (Q; V) and the fact that A is 
coercive. 

Now we only need to prove that the condition (15.41 ) is satisfied as well. 

(1) By Lemma |3~T1 we can find rj G H 1 ^; M) with divr? = \i and IMIi/i^M) — c II^IU 2 (n;V)- 

(2) Since oj + h S 2 U^ h r) G L 2 (S1;V), we can apply Lemma ED again to find t G fP(h;M) 
with divT = oj + II? h 5 2 n? h r/ and 

Ikllifi^M) - c (IMIl 2 (Q;V) + W^$,h^^-f,hV\\L 2 (Cl-V))- 

(3) Since Si is an isomorphism from i? 1 (fi;M) to M), we have p G H^fyM) with 
5i£> = r, and ||£>||.ffi(n;M) < cIMIifH^M)- 

(4) Define o- = curing g+Uj h r] G Vf+i^ 2 {Th]Y). According to Lemma 14 in ifTTl. curing o G 
Vf+iA 2 (T h ; V). So we have a G Vf+iA 2 (T h ; V). 

(5) From step (4), (16.41 ). step (1), and the fact that II 3 h is a projection, we have 

diva = divllf a t? = np h div?7 = Uf h n = //. 

(6) Also from step (4), 

-n?^ = -n r 3 h 5 2 curinJ;; e - r%s 2 n^. 

Applying, in order, Lemma [HI ( 16.311 ), ( 16.341 ), step (3), ( 16.311) , step (2), and the fact that lip h is a 
projection, we obtain 

-n? h 5 2 curini;^ = n^divSin^g = np^divn?'-^;^ 

= np ift divnj;-5x^ = n? h divn r 2 ;-r = n r 3 ft divr 
= np jh (w + np^n? ^ = w + np^n? hV . 
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Combining, we have —lip h S2<J = u. 

(7) Finally, we prove the norm bound. From the boundedness of S2 in L 2 , (16.51 ). and step (1), 

||n? h S*2np /l 77|| jL 2 (n . V ) < c||5 , 2n?^r/|| L 2( Q;V ) < c\\U^ h ri\\ L 2^ n . M ^ < c\\rj\\ H i( Q . M ) < c\\p,\\ L 2^. v y 

Combining with the bounds in step (3) and (2), this gives ||f?||i?-i(fi;M) < c (IMIl 2 (q ; v) + IMIl^^v))- 
From 4631), we then have ||curin~^£>|| L 2 (n . v) < c\\g\\ H i {n . M) < c(||w|| L 2 (Q;V) + ||/u|| L 2 (n . v) ). From 
([63]> and the bound in Step (1), ||LT? h v\\ (n-,M) < c\\rj\\ H i {n . M) < c\\n\\ L 2 (n . Y) . In view of the 
definition of a, these two last bounds imply that ||o"||l2(q ;M ) < c(||u;||£2(Q.y) + ||/u||L2m. V )), while 
||divcr|| i 2(Q. V ) = 1 1 p, 1 1 £2 (q ; v) by Step (5), and thus we have the desired bound ( 17. 11 ). ■ 

We have thus verified the stability conditions (15.31 ) and (15 -4b . and so obtain the following quasi- 
optimal error estimate. 

Theorem 34 Suppose (a,u,p) is the solution of the elasticity system ( 15.71 ) and (o~h,Uh,Ph) is tne 
solution of discrete system ( 15. 2D . where the finite element spaces satisfy the hypotheses of Theorem\33\ 
We also assume that there is r max G N such that for any h > and f € A(Th),r(f) < r max . Then 
there is a constant C, independent of h, such that 

Ik - Vh\\H(div,Q;M) + \\ u ~ '"/i||l,2(fi;V) + \\P ~ Pfrlll^V) 

<Cinf(||o--r|| ii -( diV)fl . M ) + \\u-v\\ L 2 {n . Y) + \\p - g||i2 (Q . v) ), 
where the infanum is taken over all r G V~ +1 (Th~, V), v G V~(Th', V), and q G V~(7h] V). 

8 Conclusions and future work 

In the paper, we have presented a generalization of Arnold-Falk-Winther (AFW) elements to the case 
of elements of variable order for a three dimensional domain. The proof of stability is based on the 
use of some variant of projection based interpolation operators, and a specially designed operator 
Tlf'h discussed in the text. We have proved the ^-stability for meshes with variable order under the 
assumption that there is an uniform upper bound on the highest polynomial order used. 

We plan to continue the research on several fronts. On the numerical side, we intend to implement 
and test the /ip-adaptive algorithm based on the coarse/fine grid paradigm. The code will be applied 
to a detailed study of problems with large material contrast including the streamer problem, discussed 
in lfl7l . The results obtained using the AFW elements will be compared with results obtained using 
the classical if 1 -conforming elements, in terms of memory use and CPU time. 

On the theoretical side, we will attempt to prove p-stability and, ultimately, the /ip-stability of the 
method. 
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